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Abstract: Different distortion classes of the azimuthal and cylindrical projections of the triaxial ellipsoid have been 
considered in our previous works. These projections make it possible to construct maps of the celestial bodies in plane- 
tary scale. However, for regions in the middle latitudes it is advisable to use a conic projection which was not developed 
until now. In this investigation we describe the development of three conic projections of a triaxial ellipsoid: a conic 
projection with true scale along meridians, an equal-area conic projection and a quasi-conformal conic projection. In 
derivation of the projections we use an elliptical cone tangent to a triaxial ellipsoid. The projections are calculated and 
maps in these projections are created for the first time. 
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1. Introduction 

Nowadays we have examples of regional mapping of celestial bodies in conic projections. For instance the Atlas of 
the asteroid Vesta available on the Internet (Vesta Atlas Gallery. Jet Propulsion Laboratory). Maps in this atlas were 
compiled in the Lambert conformal conic projection for the sphere despite the fact that the shape of asteroid is approx- 
imated by a triaxial ellipsoid. If the compression of shape of a celestial body is more than 40% the difference in the con- 
ic projections of sphere and ellipsoid will be significant and the reasonability of using of a triaxial ellipsoid is justified. 
The asteroid (951) Gaspra, the asteroid (243) Ida and the asteroid (433) Eros are examples of such bodies. 

For triaxial ellipsoid we assume that conic projections are projections with the meridians as a set of straight lines 
drawn from a single point, and parallels are curves constructed in accordance with the selected character of distortion. 
Herewith angles at the crossing points of meridians on the projection depend functionally on corresponding angles on 
the ellipsoid and of the cone parameters. This definition of conic projections allows us to connect the different classes 
of projections of a triaxial ellipsoid in a united system. 

In this study we present the development of three conic projections of a triaxial ellipsoid with use of an elliptical tan- 
gent cone: 

-conic projection with true scale along meridians; 

-equal-area conic projection; 

-quasi-conformal conic projection close to the conformal projection in the neighborhood of each meridian cor- 
responding to meridian section which we call the projection of the meridian section. 
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2.Derivation of conic projections 

For the derivation of conic projections we use a three-dimensional rectangular coordinate system x, y, z with the 
origin at the center of the triaxial ellipsoid and associated system of planetocentric angular coordinates 0,A on the el- 
lipsoid surface. Let us mention that planetocentric latitude is the angle between the vector r coming from the center of 
the ellipsoid and the equatorial plane and the longitude is the angle between the plane of the meridian section and the 
plane of the prime meridian The surface of a triaxial ellipsoid is described by an equation that determines the length of 
the vector r as a function of latitude and longitude for given values of the semi-axes a, b, c and calculated thereon ec- 
centricity of the ellipse of the equatorial section e x and eccentricity of the ellipse of the prime meridian section e . 

To derive formulas for conic projections we need elements of the infinitesimal trapezium on the ellipsoid corre- 
sponding to given increments of latitude dQ and longitude d/ 1 in order to compare them with the same elements on 
the projection plane (see Fig. 1,2). With an accuracy of up to infinitesimal values of a higher-order, the trapezium can 
be accepted as an infinitesimal parallelogram. 

Sides of the parallelogram are differential of the arc of the meridian and the differential of the arc of parallel. Let CO 
be the angle between them. All these parameters are expressed by Gaussian coefficients, which we obtained using the 
increment of the vector r along the curve on the surface of the ellipsoid (Nyrtsov et al., 2015). 
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Fig. 1 . An infinitesimal trapezium on the surface of ellipsoid. 

In the projection plane we use a polar coordinate system p,8 with the pole at the crossing point of the meridians 
and polar axis which corresponds to the prime meridian and the positive direction counterclockwise. Also we use a rec- 
tangular system which X pr0 j coordinate axis is directed horizontally to the right and Y pr0 j coordinate axis is directed verti- 
cally upward along the central meridian of a map with the origin at the intersection of the meridian and chosen parallel. 
Planar rectangular coordinates in projection are: 

x . = psinS , v — p — pcosS , 

pro] J pro] r' origin r ’ 

where the polar distance along the central meridian p origin = p( 0 origin ,A 0 ) is equal to distance between an origin of 
polar coordinate system and the origin of planar rectangular coordinate system. 
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Fig. 2. Infinitesimal trapezoid (parallelogram) in conic projection. 


For the derivation of conic projections we obtain the tangent cone parameters and relation of the angle 8 at the vertex 
of the cone and longitude on the ellipsoid. The formula for calculation of the angle 8 : 

4 

8 = j* OcdX , where OC is the function of longitude and the cone parameter. 

A 0 

We derive formulas of the polar distance for all three offered projections. 

For conic projections with true scale along meridians 

p t = \4EdO + L 

0j 


For equal-area conic projection 


Pi = 


For the conic projection of the meridian section: 


p i = L exp 


L 2 +- j^EGsincod# 

4e 


a J 


dO 


y sin co j 


Here L is the length of the tangent line segment from the tangent point to the vertex of the cone; O tg - the latitude of 
the tangent point. 

After all we calculate distortions in conic projections. Ellipses of distortion (Tissot indicatrix) for conic projections 
are shown in Figure 3. Meridians and a line of a tangency of a cone and triaxial ellipsoid do not depend on character of 
distortions unlike to parallels. Ellipses are connected to the cartographic grid with the parameters of the asteroid 433 
Eros and with the direction of counting of longitude goes to the west. 
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Fig. 3. Tissot indicatrix for a) conic projection with true scale along meridians, b) equal-area conic projection, c) conic projection of the me- 
ridian section. Red color line is a line of tangency of a cone. Blue color trapeze is an area of map depicted in the Fig. 4 


3. Mapping of asteroid 433 Eros in conic projections 


The asteroid 433 Eros was taken for approbation of conic projections of the triaxial ellipsoid. It is approximated by a 
triaxial ellipsoid with semi-axes a = 17 km, b = 5.5 km, c = 5.5 km (Report of the IAU Working Group, 2009). Maps 
were created in conic projection with true scale along meridians, equal-area conic projection, conic projection of the 
meridian section for region with craters Casanova and Pelleas from 40° to 70° latitude and from 210° to 250° longitude. 
Central meridian of map is 229°, latitude of center of a map - 55°. Maps are shown in Fig. 7. To show the asteroid’s sur- 
face we used photo mosaic prepared by Philip J. Stooke and available through the Small Bodies Node of NASA’s Plan- 
etary Data System. 
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Fig. 4. Maps of asteroid 433 Eros in a) conic projection with true scale along meridians, b) equal-area conic projection, c) conic projection of 
the meridian section. 


4. Conclusion 


In the article we described derivation of formulas for three conic projections of triaxial ellipsoid for mapping of the 
territories in the middle latitudes of celestial body. For the triaxial ellipsoid as a surface which can be projected on a 
plane without distortions we use a direct elliptic cone tangent to ellipsoid. 

The polar coordinates system is used for derivation of projections. The polar angle is obtained by integration on a lon- 
gitude of an increment of this angle depending on an increment of a longitude. The polar distance for all three offered 
projections is obtained by integration on a meridian that is with the account of an increment of latitude at the fixed value 
of a longitude. The calculation of the coordinates of the nodes of the cartographic graticule and values of distortions of 
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offered projections are realized in the GIS Research Centre Site by IG RAS (GIS Research Centre Site, Conic projec- 
tions, 2017). 
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